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Abstract. We review cosmological backgrounds of gravitational waves with a partic- 
ular attention to the scientific potential of the eLISA/NGO mission. After an overview 
of cosmological backgrounds and detectors, we consider different cosmological sources 
that could lead to an observable signal. We then study the backgrounds produced by 
first-order phase transitions and networks of cosmic strings, assessing the prospects for 
their detection. 
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1. Introduction 

Because of the weakness of the gravitational interaction, gravitational waves (GW) 
propagate freely during essentially all the cosmic history. Once produced in the early 
universe, they simply red-shift with the cosmological expansion, retaining their ini- 
tial frequency dependence. This allows GW detectors to directly probe scales, epoch s 
and physical phenomena tha t are n o t access i ble in any other way (Refs. lAlienl (l996), 
iMaggiord d2000h . iBuonannol d2003h. iHoead (l2006h). 

In this contribution, based on Bi netruy et all (120121) . we review cosmological GW 
backgrounds with a particular atte ntion to those backgrounds t hat could be observ- 
able by the eLISA mission (Refs. |X maro-Seoane et all (l2012ah . lAmaro-Seoane et all 



d2012bh ). eLISA (also named NGO by ESA) is the new, European version of LISA that 



is now being proposed Q. We start with an overview of cosmological backgrounds and 
detectors in Section [2] In Section [3j we discuss several cosmological sources of GW. 
We then study the backgrounds produced by first-order phase transitions (Section @]) 
and cosmic strings (Section©, assessing the prospects for their detection. We conclude 
in Section [6l 

Notation: we work in units where c — ks — 1 • 



2. Cosmological Backgrounds and Detectors 

A cosmological source of GW typically consists of many uncorrected and unresolved 
events, which produce a stochastic GW background. In general, such a backg round 
is expected to be nearly Gaussian, isotropic, stationary and unpolarized (Refs. lAllen 



'See http://www.elisa-ngo.org/ for more details. 
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( 19961) . iMaggiord (|200Q|) ). Its main property is then its frequency spectrum. The quan- 
tity that is usually considered to characterize cosmological backgrounds is the spectrum 
of GW energy density per logarithmic frequency interval, normalized to the critical en- 
ergy density today 

h 2 dp 



h^a gw {f) = --7TZT7 W 



<„2 r> <f\-'±_ " H 8 W 

Pc dlogf 

where / is the present-day GW frequency, p gw the present-day GW energy density, 
p c - 3H 2 /(SnG) is the critical energy density today, and h - 0.72 parametrizes the 
small uncertainty in the value of the Hubble constant today (H = 100/jkm/s/Mpc). 

There are currently a number of upper bounds on cosmological GW backgrounds, 
in very different frequency ranges. The amount of GW that is present at the time of Big 
Bang Nucleosynthesis (BBN) is con strained not to s poil the successful predictions for 
the light element abundances (Ref. Maggiore! (120001) ). Although this bound applies to 
the total GW energy density p gw , neglecting the unlikely possibility of a GW spectrum 
with a very narrow peak one can take as a rule of thumb h 2 D. gw < 10~ 5 for / > 10 -10 Hz 
(/ ~ 10" 10 Hz corresponds to the comoving Hubble scale at the time of BBN). Simi- 
larly, the total GW energy density that is present at the time of decoupling of the Cosmic 
Microwave Background (CMB) is constrained by precise measurements of the CMB 
power spectrum (Ref. Smit h et all (l2006bl) ). The current bound is also of the order of 
h 2 Q. gw < 10" 5 ,butit covers a larger frequency range, / > 10 15 Hz. At the lowest ob- 
servable frequencies, / ~ 10~ 18 - 10~ 16 Hz, cosmological GW are constrain ed by their 
contrib ution to the large-angle fluctuations of the CMB temperature (Ref. Maggiore 
(2000)). The strongest bound is obtained at / ~ 10" 16 Hz: h 2 Q. gw < few x 10 -14 . The 
sensitivity to cosmological GW in this frequency range ma y improve significantly in th e 
future with measurements o f the CMB polarization (Refs. lKamionkowski et al.l(ll997h . 
Seljak & Zaldarriaga] dl997l) ). GW back grounds are also constraine d by the very ac- 



curate timing of msec pulsar arrays: Ref. Ivan Haasteren et all (120111) recently obtained 
the upper bound h 2 Q. gw < 6 x 10~ 9 at / = 4 x 10~ 9 Hz (up to a weak dependence 
on the slope of the spectrum). The sensitivity of pulsar timing arrays t o stochastic 
GW b ackgrounds is expected to improve significantly in the future (Ref. Jenetetal 



(2006)). Finally, the S5 run of LIGO allowed to put the upper bound (Ref . lAbbott et al 
hm% ) h 2 Q gw < 3.6 x 10~ 6 at / « 100 Hz, which improves upon the BBN and CMB 
bounds in this frequency range. Again, the sensitivity to stochastic GW backgrounds 
in this frequency range will improve significantly in the near future with the advanced 
ground-based interferometers. 

For ground-based interferometers, the best way to disentangle a stochastic back- 
ground fr om the instr u mental noi s e is to cross-correlate the output of different detec- 
tors (see lAllenl ( 1996 ). Maggiore (2000) and references therein). This is not possible 
for eLISA since there is a single instrument. In this case, the spectral density of the 
signal must be larger than the one of the noise. However, the fact that eLISA is de- 
signed to probe lower frequencies (/ ~ 10~ 4 - 10~ ! Hz) is an important advantage, 
becau se the spectral de nsity of the signal Sh(f) is related to Q. gw as D. gw oc f 3 Sh(f)/H 2 
(Ref. iMaggiorel (|2000)). Thus, for a given sensitivity to S /,(/), the sensitivity to the 
GW energy density quickly improves as the frequency decreases. Compared to the 
original LISA mission, it is more difficult for eLISA to disentangle a stochastic back- 
ground from the instrumental noise, all the more because the eLISA Mother-Daughter 
configuration, providing only tw o measurement arms, does not allow to use Sagnac cal- 
ibration (Refs. lTinto et al.l ( 2001 ).lH o gan & Bender] ( 200lh ). Luckily as we will see, for 
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sources like phase transitions and cosmic strings, the spectral dependence of the signal 
is well predicted and may allow to distinguish cosmological backgrounds as long as 
they lie above the sensitivity curve. In Sections @] and [5] we consider the sensitivity 
curve of eLISA computed from the expected instrumental noise an d after removing the 
confusion noise generated by unresolved galactic binaries (Ref. Amaro-Seoane et al. 
(l2012al) ). 



3. Cosmological Sources of GW 

A cosmological source emits GW with a characteristic wavelength A* that is smaller 
than the Hubble radius H~ Y at that time: A* = e, H~ l with e* < 1. This wavelength is 
then stretched by the expansion of the universe from the time of production until now. 
For GW produced during the radiation era when the plasma temperature is T t , and 
assuming a standard adiabatic thermal history for the evolution of the universe after 
GW production, the characteristic GW frequency today can then be written as 

^ 1.6x10-* Hz |M« 

where g t is the number of relativistic degrees of freedom at the temperature r*. The 
parameter e* < 1 depends on the dynamics of the particular GW source under consid- 
eration. For a first-order phase transition, one may have for instance e* ~ 10~ 3 - 1. In 
that case, Eq. © shows that GW produced around the electroweak scale may fall right 
into the mHz frequency band of eLISA. In the case of cosmic strings, GW are produced 
continually during the cosmological evolution (over a wide range of values of T t ), so 
that the present-day spectrum covers a very wide range of frequencies. We will discuss 
these two sources in Sections |4] and [5] 



Several other cosmological sources of GW are reviewed in detail in Bin etruy et al. 

One source that is very well motivated by other observations is inflation. 
During in flation, quantum fl uctuations of the graviton field are parametrically ampli- 
fied (Ref. iGrishchukl (Il975h ) into tensor perturbat ions at supe r -Hubb le scales by the 
quasi-exponential expansion of the universe (Ref. IStarobinskvl dl979h \ These tensor 
perturbations become standard GW when they re-enter the Hubble radius in the course 
of the post-inflationary evolution, leading to a very broad spectrum today. If inflation 
occurs at sufficiently high energy scales (not much below Ej n fi ~ 10 16 GeV), the result- 
ing GW in the frequency range / ~ 10~ 18 - 10" 16 Hz ma y be indirectly detected in th e 
future through their effect on the CMB polarization (Refs. iKamionkowski et all dl997D . 
Seljak & Zaldarriaga| (Il997l) ). In the vast majority of cases, the current CMB bound on 



the inflationary GW leads already to the upper bound h 2 Q gw < 10 15 on their amplitud e 



at frequencies accessible by direct detection experiments, see e.g. ISmith et all ((2006a). 
This is much below the reach of eLISA and ground-based experiments. 

Under certain circumstances, however, the amplitude of the inflationary GW back- 
ground may be enhanced in the direct detection frequency range, while satisfying the 
CMB constraint at low frequencies. This occurs in some inflationary models where par- 



ticles a re abundantly produced during inflation, see lCook & Sorbol ( 2012 ). lBarnaby et all 



(2012). Another possibility is that the early universe was dominated during some epoch 
after inflation by a fluid stiffer than radiation, i.e. with an equation of state w > 1/3 (see 
iBuonannol d2003h and references therein). A GW background with an enhanced am- 
plitude at high frequencies may also be produced in some alternatives to inflation like 
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the pre-big-bang scenario (see lBuonannol d2003l) and references therein). The GW pro- 
duced by these processes may be observable by direct detection experiments, preferen- 
tially ground-based experiments since they probe higher frequencies. 

Another privileged epoch for GW production is the end of inflation, when the 
potential energy density driving inflation is converted into the thermal bath of the 
Hot Big Bang in the course of reheating. In many inflationary models, reheating 
starts with an explosive and non-perturbative decay of the inflaton condensate into 
very large fluctuations of itself and other bosonic fields coupled to it - a process called 
preheating. The large field flu ctuations amplified by prehea t ing can produc e a sig - 
nificant amount of G W (Refs. |K hlebnik ov & Tkachevl I dl997h. lEasther & Liml (2006), 
Easther et al. (2007), Garcia-Bellido & Figueroa (2007), Garcia-Bellido et al. (2008), 



Dufaux et all (120071) . iDufaux et all (12009T) . Dufaux et all (I2010h . iBrax et all (1201 lh ). 
The peak frequency of these GW today depends on the energy scale of inflation: it 
may fall into the frequency range of ground-based interferometers if E^fi < 10 11 GeV 
and in the eLISA frequency range if < 10 7 GeV. Thus the GW produced by pre- 
heating may be observable if inflation occurs at a sufficiently small energy scale, which 
is complementary to the GW produced by inflation itself. GW can also be produced 
by the non-perturbative decay of scalar fi eld condensates other than the inflaton, in par- 
ticular in supersymmetric theories (Ref. IDufauxl (12009b ). In this case, the peak GW 
frequency today is not directly related to the energy scale of inflation but rather to the 
scale of supersymmetry breaking. 

Several potential sources of GW are related to the spontaneous breaking of sym- 
metries during early universe phase transitions. Besides first-order phase transitio ns and 



cosmi c strings, another possibi l ity is provided by unsta ble domain walls (see e.g. lGleiser & Roberts 
(ll998h . lTakahashi et all d2008h . iHiramatsu et al.l(l2010l) ). These could produce GW po- 
tentially observable by eLISA, depending on the details of the underlying high-energy 
physics model. GW can also b e produced by the self-ordering o f a scalar field after a 
global phase transition (see e.g. lKraussldl992h . lFenu et aT1 (2009). K rauss et~aT1 ( 201C)l) ). 
An observable amount of GW for eLISA would require a scalar field with a very large 
expectation value. 

As another example, GW may also be produced by l arge density perturbations 
leading to the formation of primordial black holes (PBH), see lSaito & Yokoyamal (t2010h 
and references therein. The resulting spectrum is peaked at a characteristic frequency 
that depends on the PBH mass. For instance, eLISA could probe the formation of PBH 
in the mass range 10 22 - 10 25 g, which could be dark matter candidates. 

In general, different cosmological sources lead to GW backgrounds with differ- 
ent frequency dependences, which may allow to distinguish them from each other and 
from astrophysical and instrumental backgrounds. It is thus important to determine the 
frequency dependence of the signal as accurately as possible. 



4. First-Order Phase Transitions 

In the course of its adiabatic expansion, the universe may have undergone several phase 
transitions driven by the temperature decrease. They occur because the effective poten- 
tial of scalar fields (e.g. the Higgs field) depends on their interactions with other parti- 
cles in the thermal plasma. The nature of the phase transition depends on the underlying 
high-energy physics model. If it is first-order, then at a certain critical temperature the 
scalar field becomes trapped in a false vacuum state separated from the true vacuum by 
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a potential barrier. In that case the phase transition proceeds by tunneling across the 
potential barrier, leading to the nucleation of true vacuum bubbles which quickly grow 
into the false vacuum phase and collide. This is a violent and highly inhomogeneous 
process that can emit a significan t amou nt of GW, as first discussed in Witten (1984), 
Hosanl (l!986h . lTurner & Wilczekl (Il990h . 



For phase transitions that occur at the electroweak scale, we saw below Eq. © that 
these GW may fall right into the eLISA frequency range. We will therefore focus on 
these in the following. In the Standard Model of particle physics, the electroweak phase 
transition is not first-order and it is not expected to lead to any significant production 
of GW. However, it can be first-order in well-motivated exten sions of the Standard 
Model, in particular those addressing the hierarchy problem (see Binet ruy et all (1201 2|) 
for references to specific models). An observable GW signal from an electroweak phase 
transition would thus provide a clear indication of physics beyond the Standard Model. 

GW from first-order phase transit i ons are produced by two ma i n processes: bubble 
collisi o n (Refs.lKosowskv et al.ldl992h.lKosowskv & Turned (1 1993hJKamionkowski et al. 



( 1994), ICaprini et all (120081) . iHuber & Ko nstandin ( 2008a), Caorini et all (l2009al)) and 



magneto-hydrodynamic turbulence (Refs. Kamionkowski et al. ( 1994), Kosowsky et al. 



(120021) . iDolgov et al.ld2002l) . ICaprini & Purred (l2006lklCaprini et al.l(l2009bl)). T he cal- 
culation of the resulting GW background is reviewed in detail in lBinetruy et all (120121) . 
It depends on five main parameters: the temperature of the phase transition T*, its du- 
ration p~ l {fl is related to the rate of bubble nucleation and should be larger than the 
Hubble rate, /? > //*), its strength a - p vac /p r ad (Pvac i s the vacuum energy density and 
p ra d the radiation energy density at the time of the phase transition), the velocity of the 
bubble walls vj,, and an efficiency factor k that measures the fraction of vacuum energy 
that is converted into gradient energy of the scalar field and bulk kinetic energy of the 
fluid (as opposed to thermal energy). These five parameters are not independent of each 
other and must be computed in the context of a specific high-energy physics model: T*, 
P and a depend on the scalar field potential at finite temperature, while v# and k depend 
on the propa gation of t he bub ble wall in the surrounding plasma. 

In lBinetruy et al.1 (l2012h . we calculated the GW background from first-order phase 
transitions by combini ng consistently the latest res ults available for GW production by 



bubble collision (Ref. Huber & Konstandin (2008a)) and turbulence (Ref. Caprin i et al 



(|2009bh ). We also used the model of bubble propagation recently developed by lEspinosa et al 



(120101) to evaluate the parameters Vb and k (these parameters were previously evalu- 
ated under the particular assumption of Jouguet detonation, which is not justified in 
most cases). This ca uses the GW signal to depend on the parameter n introduced in 
Espinosa et all (120101) . which characterizes the friction exerted on the bubble walls by 



the surrounding plasma. 

We show the resulting GW backgrounds in t wo specific high - energy models with 
a first-order phase transition in Figs. [Tl and |21 (see iBinetruy et all d2012l) for a model- 
independent discussion of the detection prospects). The frequency dependence of the 
signal is well determined in each case. It varies as Cl gw oc f 3 at low frequencies, while 
the high-frequency tail is given by the combination of Q. gw oc f~ x from bubble collision 
and Q. gw oc / _5 ^ 3 from turbulence. 

In Fig. Q]we show the GW spectra f or different values of the parameters in a model 
studied in Huber & Konstandin (2008b) where the Higgs potential includes a dimen- 
sion six operator. For small values of a (lower spectra), one clearly distinguishes two 
different peaks, a low-frequency one due to bubble collision and a high-frequency one 
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f (Hz) 

Figure 1. GW spectra from the electroweak phase transitio n in a model where 
the Hi ggs potential includes a dimension six operator (Ref. Huber & Kon standinl 
(2008b)), compared to the eLISA sensitivity curve. The GW spectra correspond to 
different values of the model parameters with a varying from 0. 128 (lower spectrum) 
to 2.268 (upper spectrum) and r\ — 1 . 



due to turbulence. The relative contribution of these two sources depends on a and 77. 
The GW amplitude increases for values of the model parameters leading to larger val- 
ues of a, but the peak of the spectrum is then also shifted towards smaller frequencies. 
Furthermore, the phase transition cannot be arbitrarily strong in this model (otherwise 
the transition never completes). As a consequence, the amplitude of the GW signal 
saturates below the eLISA noise curve. 

More promising is the phas e transition of the ra dion field i n the Randall-Sundrum 
model, which has been studied in Randall & Servant! ( 2007 ) and Konstandin et all ( 20101) . 
In that case, the phase transition is very strong, a » 1 . The GW signal is then dom- 
inated by bubble collision with k - 1 and Vb — 1. Furthermore, the amplitude of the 
GW background becomes independent of a for a » 1 . The signal then depends only 
on two parameters: /?///* and T t , which can take a wide range of values in this model. 
In Fig. f2J we show the GW spectra for T t = 100 GeV (left spectra) and T t = 10 4 
GeV (right spectra). In both cases, the GW amplitude increases with the duration of 
the phase transition . A large signal-to-noise ratio for eLISA is obtained for slow 
enough phase transitions with T t > 10 3 GeV. For such large values of T*, the model is 
beyond the reach of the LHC, so eLISA might be the only way to probe it. 



5. Cosmic Strings 

Cosmic strings (Ref. Vilen kin & Shellardl ( 1994 )) are linear concentrations of energy 
with a cosmological size that arise in a variety of high-energy physics mode ls, both in 
field th eory (grand unification, supers ymmetry, ...) and in string theory (Ref. iPolchinski 
(2004)). They are relics from phase transitions occurring at the end of inflation or dur- 
ing the thermal evolution of the universe. Once produced, a network of stable cosmic 
strings evolves towards a self-similar regime characterized by a continuous energy loss: 
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f (Hz) 

Figure 2. GW spectra from the radi o n phase transition in the Randall-Sundrum 
model (Refs. iRandall & Servant! (120071) . iKonstandin et al.l (120101) ). compared to the 
eLISA sensitivity curve. The temperature of the phase transition is T* - 100 GeV 
for the left/plain spectra and T* = 10 4 GeV for the right/dashed spectra. In both 
cases, /3/H t = 6 for the upper spectrum and 15 for the lower one. 



when long string segments cross each other and reconnect, they form smaller cosmic 
string loops, which oscillate relativistically and decay away by emitting gravitational 
waves. The net result is that a significant fraction of the cosmic string energy density 
is continually converted into GW via the production of loops, from t he very early uni - 

verse until the presen t epoch. The resulting G W ba ckground (Refs. Vilerikinl (11981 ), 

Hogan & Rees ( 1984), Vachas pati & Vilenkid(fl985l) . ICaldwell & Alien! dl992n . lSiemens et al 
d2007l) . lDePies & Hoganl d2007l) ) covers a very wide frequency range and can be looked 
for by different experiments, in particular pulsar timing arrays, ground-based inter- 
ferometers and eLISA. In addition to this background, the GW signal from cosmic 
strings includes also strong and infr equent bursts produc e d by s pecific loop configu- 
ration s called cusps and kinks (Refs. iDamour & Vilenkir] (l200lh . iDamour" & Vilenkin 



(2005)). These GW bursts can be looked for individually and should not be included in 
the computation of the stationary and nearly Gaussian background. 

The GW signal from cosmic strings varies with two fundamental parameters that 
depend on the underlying high-energy physics model: the string tension p (which is 
equal to their energy per unit length in the simplest models that we consider below) and 
the probability p that two string segments reconnect when they cross each other. The 
dimensionless parameter Gp (where G is Newton constant) must be sufficiently small 
to satisfy observational constraints on the gravitational effects of cosmic strings (CMB, 
gravitational lensing, ...)q The reconnection probability is usually p = 1 for cosmic 
strings in field theories, but it can be smaller for cosmic strings in string theory, which 
incr eases the string energy den sity p s . Simple arguments suggest that p s oc l/p (see 
e.g. IDamour & Vilenkin (2005) ), but a mo re complicated dependence has also been 
found in e.g. Avgoustidis & She llardl d2006l) . We consider for definiteness p s oc l/p be- 



2 The production of particles by cosmic strings m ay also lead to complementary constraints on the string 
tension, see e.g. iDufauxl j2012al) . lDufauxl j2012bh . 
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low, but the results can be directly generalized to the case where p s x 1 / f(p) with an ar- 



bitrary function f(p) by simply replacing p by f(p) in Figs. [3] and HI see lBinetruy et al. 
(2012) for more details. 

The GW background from cosmic strings depends also strongly on the typical 
size of the loops when they are produced, which is still uncert a in by several orders 
of magnitude. The latest simulations (Ref. iBlanco-Pillado et all (1201 lb ) indicate that 
the loops are produced with a size about ten times smaller than the Hubble radius. 
The loops are then long-lived compared to the Hubble time. F or definiteness, we wi ll 
focus on this case below ; other cases are discussed in detail in Binetru y~et all (120121) . 
Also uncertain is the typical GW spectrum emitted by each individual loops, which 
depends in particular on the presence of cusps and kinks. However, we showed in 
Binetruy et al.1 (1201 2|) that this uncertainty does not strongly affect the predictions for 
the GW background once the total power (P gw ~ 50G/J 2 ) emitte d in GW by each loo p is 
specified. Compared to previous works, we also improved in lBinetruy et al.l d2012l) the 
modeling of the cosmological evolution, because it affects the frequency dependence 
of the GW background. 
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Figure 3. GW background from cosmic strings compared to observational sensi- 
tivities, for different values of the cosmic string parameters in the case of long-lived 
loops (see the main text for details). 



The resulting GW background for long-lived loops is shown in Fig. [3] for differ- 
ent values of Gp and p. The frequency dependence of the signal is well determined. 
The high-frequency part of the spectrum (which is produced during the radiation era) 
is almost flat over many decades of frequency (small deviations from an exactly flat 
spectrum arise because of the variation of the number of relativistic species during the 
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expansion of the universe). The spectrum slightly increases at lower frequencies and 
then falls off as Q. gw oc f 3 / 2 at very small frequencies. 




Figure 4. Regions of the parameter space of cosmic sttings that can be probed by 
different observations (at the right of each curve), in the case of long-lived loops (see 
the main text for details). 

An important feature of the GW background from cosmic strings is that, because 
it covers a very wide frequency range, it can be observable by different experiments. 
We show the sensitivity of eLISA, ground-based interferometers and pulsar timing ar- 
rays in Fig. [3] For ground-based interferometers, we consider the upper bound set by 
the LIGQ S5 run and the expected sensitivity of Advanced LIGO, both taken from 
Abbott et~ai1 (I2009I) . For p ulsar timing experiments, we consider the upper bound ob- 
tained in ljenet et all (|2006T ) and the future sensitivity of the complete Parks Pulsar Tim- 
ing Array (PPTA) estimated in that reference. 

In Fig. HI we show the regions of the parameter space that can be probed by these 
experiments, still in the case of long-lived loops. The darkest region is already excluded 
by current pulsar bounds. In a large part of the parameter space, the GW signal could 
be simultaneously observed by eLISA and at least one other experiment. This would 
certainly improve the characterization of the signal and the ability to distinguish it from 
other stochastic backgrounds. Other regions of the parameter space are only accessible 
to eLISA, which could thus open a new window on the cosmic string parameters. In 
particular, eLISA would be able to probe most of the parameter space for cosmic strings 
that arise in the simplest string theory models. 

6. Conclusion 

We have investigated the scientific potential of the eLISA mission regarding the de- 
tection of cosmological GW backgrounds. We have considered different cosmological 
sources that could lead to an observable signal, with a special emphasis on first-order 
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phase transitions and cosmic strings. Each of these two sources produces a GW back- 
ground with a distinct and well determined frequency dependence, which allows to 
distinguish the signal from other cosmological or astrophysical backgrounds and from 
the instrumental noise. In both cases, the signal may lie well above the eLISA noise 
curve, depending on the microphysics of the source. The detection of a cosmological 
GW background by eLISA is of course not guaranteed: it could just constrain cosmo- 
logical sources. However, such a detection is at least possible, and this would be a 
major discovery for our understanding of the very early universe and the underlying 
high-energy physics. 

Acknowledgments. It is a pleasure to thank Pierre Binetruy, Alejandro Bohe and 
Chiara Caprini for collaborations and the organizers of the ninth LISA symposium for 
a very interesting conference. 
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